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Abstract 

We present a complete calculation of the 0{a s ) QCD corrections to the 
width of the decay H + — * tb within the Minimal Supersymmetric Standard 
Model. We find that the QCD corrections are quite important, and that the 
supersymmetric QCD corrections (due to gluino, t and b exchange) can be 
comparable to or even larger than the standard QCD corrections in a large 
region of the supersymmetric parameter space. This is mainly due to the effect 
of large left -right mixings of stop (t) and sbottom (6). This could significantly 
affect the phenomenology of the H + search. 
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1 Introduction 



The existence of a charged Higgs boson H + would be a clear indication that the Stan- 
dard Model must be extended. For example, the Minimal Supersymmetric Standard 
Model (MSSM) Jl| with two Higgs doublets predicts the existence of five physical 
Higgs bosons h°,H°,A°, and H ± || |3|]. If all supersymmetric (SUSY) particles are 
heavy enough, H + decays dominantly into tb above the tb threshold 0, f|]. In refs. |5|. |6[ 
all decay modes of H + including the SUSY-particle modes were studied in the case 
that the SUSY-particles are relatively light: it was found that the tb mode remains 
important even in this case (though the tb mode can be dominant in a wide range 
of the MSSM parameters). Thus it is important to calculate the QCD corrections to 
the tb mode as they could significantly affect the phenomenology of the H + search. 
The standard QCD corrections to the tb mode were already calculated 0: they can 
be large (+10% to -50%). There also exist calculations of the SUSY-QCD correc- 
tions within the MSSM || |§. However, in ref. || the squark-mixing was neglected. 
The calculation in ref. is incomplete as the wave function and mass renormaliza- 
tions were omitted. A calculation of SUSY-QCD corrections to the related process 
t — > H + b was done recently in ref. []T0| . 

In this paper we present a complete calculation of the 0(a s ) QCD corrections to the 
width of H + — > tb within the MSSM. We include the left-right mixings of both the 
ti,R squarks and the bL,R squarks. We adopt the on-shell renormalization scheme. 

2 QCD one— loop contributions 

The one-loop corrected amplitude of the decay H + (p) — > t{k t )b{ki) (p = k t + ki) can 
be written as 

M = ii(Y l P R + Y 2 P L )b (1) 
with Pr,l = |(1 ± 75) and the one-loop corrected couplings: 

Y i = y i + 5Yt 9) +SYf 9) (i = l,2), (2) 
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where i/i are the tree-level couplings corresponding to Fig. la: 

Hi = —j^- — m b tan (3 = h b sin j3 , 
y2m w 

y 2 = —tM- — m t cot (3 = h t cos(3, (3) 
y2m w 

with g being the SU(2) coupling. 8Y^ and 5Y® are the contributions from gluon 
and gluino exchanges, respectively (as shown in Figs, lb and lc). 
The tree-level decay width is given by: 

T tree {H+ _> fb) = _^[« + - m \ - m 2 b \y\ + y 2 2 ) - Am t m byi y 2 ] , (4) 
Ibirm H+ 

where k = n(mj I+ , m 2 , ml), k(x, y, z) = ((x — y — z) 2 — Ayz) 1 / 2 , and N c = 3. 



The vertex corrections due to gluon and gluino exchanges at the vertex (Fig. lb), 
SYf ' 9 ^ and 5Y^ v ' g \ respectively, are given by: 

5(Y 1 P R + Y 2 P L )^ = ^{2[B (mlO,m 2 t ) + B (mlO,m 2 b )-r 

-{m 2 H + -m 2 -ml)C (\ 2 ,rn 2 ,m 2 b )](y 1 P R + y 2 P L ) 
-2miC 1 (A 2 ,m^,m^)[(m t ?/ 1 +m b y 2 )P R + (m t y 2 + m 6 yi)P L ] 
-2m b C 2 (\ 2 ,m 2 ,m 2 b )[(m t y 2 + tu^Pr + {m t y x + m b y 2 )P L ] } , 

5(Y 1 P R + Y 2 P L f^ = ^ { 2Gij [-m- 9 C (m 2 ,ml,ml){(a LR ) t ,P R + (a^-i^} 



47T 

+m t C 1 (mj,m~.,m 2 .){(a LL )i j P L + {a RR )ijP R } 
+m b C 2 (m 2 § , m 2 .,m 2 bj ){(a LL ) ij P R + {a RR )ijP L ]\ } . (5) 
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with Cf = 4/3 and 



/ cos &i cos 9 b — cos $i sin ^ \ / — cos 9~ t sin 9~ b — cos % cos 9 b \ 

aLL ~ \ -sin0 t -cos0 6 sin t - sin ^ J' aLfl ~ \ sin 9~ t sin 6 sin0 t ~cos0 s J ' 

/ — sin 9i cos 9 b sin 0^ sin 9 b \ ( sin 0^ sin 9 b sin 0^ cos 9 b \ , , 

aRL = ^ -cos0 t ~cos0 g cos%sin0g J ' ^ = ^ cos t ~ sin S cos0 t ~cos0 s J ' W 

Gjj are the tree-level couplings of H + to = 1, 2) reading: 

^ g ff t / m 2 tan /3 + m 2 cot (3 - m 2 ^ sin 2/3 m b (A b tan(3 + fi) \ , b ^ 

13 V2m w { m t (A t cot(3 + ii) 2m t m b / sin 2/3 J l ^ j " 

(7) 



Here R q (g = t or b) is the g-mixing matrix 



Rt 



cos 9g sin 6q 

— sin 6q COS 6q 



1,2; a = L,i?) 



relating the squark states g^ and q R to the mass-eigenstates q~\ and g 2 < m q 2 ) '■ 
q"i = R\Jla- R? diagonalizes the squark mass matrix H]: 



m 
rn 



LL 
2 

RL 



m 
m 



LR 
2 

RR 



(R 



m| 











R q , 



rn 



'12 



where 



rn 



LL 



Ml + m 2 q + m| cos 2(3(I q - Q 9 sin 2 6 W ) , 



m 



RR 



M 



+ m 2 + m| cos 2/?Q g sin 2 



n? 



LR 



m 



RL 



{u 

m t {A t - //cot/3) (<? = £) 
mb(y4b — /x tan /?) (g = b) 



(9) 

(10) 
(11) 
(12) 



As usually, we introduce a gluon mass A for the regularization of the infrared diver 

nd C 
[1,?" 



gence. Here we define the functions B , Bi, C , Cx, and C 2 as in [11], [12 

d D g 



[5 ,rai](fc 2 ,m 2 ,m ) = /Hr71 2 ' /u2 2T ( 13 ) 

j m 2 (q 2 — mo){{q + k) 2 — m 2 ) 

Co,kfCi — k^C2 (t?t.q, 77T. 2 , w 2 ) = / — ^--r^ 2\77 , ; \! ~2\77 TT5 2\ • 
6 J J 27H (g 2 - mg)((g + k t ) 2 - mf)((g - fcg) 2 - m^) 

Here fc^ and ^ are the external momenta of t and 6, respectively. The parameter r in 

eq. (0) and following equations shows the dependence on the regularization: r = 1 for 

dimensional regularization and r = for the dimensional reduction (DR) fI5| . The 

dependence on r, however, disappears in our final result. 



Now we turn to the quark wave-function renormalization due to the graphs of 
Fig. lc. The two-point vertex function for gg can be written as: 

ff(l + K(k 2 )P L + U R (k 2 )P R ) - K + mk 2 )P L + mk 2 )Pn). (14) 

Here we have S|(fc 2 ) = Y, q R (k 2 ) = Y, q (k 2 ). The correction to the amplitude from the 

wave-function renormalization has the form: 

S(XiPr + Y 2 P L ) {W) = -|(nl(m 2 ) + n R (m 2 )) yi P R - i(n^(m 2 ) + Tl L (m 2 ))y 2 P L 
+ {m t t t {m 2 ) - m 2 n*(m 2 ) + m b ± b (m 2 ) - m 2 b U b (m 2 ))( yi P R + y 2 P L ) , 

(15) 
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with IF = ±(rif, + iifj) and X ee ||. The explicit calculation yields: 
(nl(fc 2 )P L + U%(k 2 )P R )^ = ^[-2B x {k\mlX)-r], (16) 



(Ul(k 2 )P L + U R (k 2 )P R )^ = -^[2(cos 2 9^P L + sm 2 e^P R )B 1 (k 2 ,mlml 

+ 2(sin 2 9 q P L + cos 2 Q q P R )B x {k 2 , m 2 , m| 



s9(A .2 )(g) = ^^[^(fca,^, A 2 )-2r], (17) 
2,(^)0) = ^[ m -sin2^(5 (A ; 2 ,m 2 ,m 2 i )- J B (A; 2 ,m 2 ,m 2 2 ))]. (18) 

Finally, there are additional corrections 5Y^ by the renormalization of the quark 
masses in the couplings of eq. (|3|) (In the DR scheme equivalent corrections are nec- 
essary as one takes the physical masses of the quarks as input): 

SY^ = 5 yi = —f^- — 5m b tan (5 , 
V2m w 

5Y^ = 5y 2 = -J— 5m t cot(3, (19) 
V 2m w 



with Sm q = 5m q 9 ^ + 8m^' , 

H 9) = -^[2m q (B (mlO,ml)-B 1 (mlO,m 2 q )- r -)}, and 

Sm q g) = s — - [ sm2e q mg(B (m 2 m 2 q ,m 2 qi ) - B (m' t ?,m|,m? )) 



+ m (? ( J B 1 (m 2 ,m 2 ,m 2 i ) + J B 1 (m 2 ,m 2 ,m 2 2 ))] .(20) 

Taking all contributions eqs. ([| [|, [15|, [19]) together, we get the one-loop corrected 
couplings Yi = yi + 5^ = ^ + SY^ + + 8Y} w) with contributions to tf^W'W-M 
from gluon and gluino exchanges (see eq. @)- It can be readily seen that they are 
ultraviolet finite but still infrared divergent. The one-loop corrected decay width to 
0(a s ) is then given by 

v(H + ^tb) = ^ [« - m " - m l) (yl + vl + 2yiRe(gyo + 2y 2 Re(^r 2 ; 



-4m f m fe ( yi j/ 2 + 2/^(5^) + l/ 2 Re(5F 1 ))] . (21) 
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3 Inclusion of the gluon emission 



For the cancellation of the infrared divergencies (A — ► 0) it is necessary to include the 
0(a s ) contribution from real gluon emission as shown in Fig. Id. 
The decay width of H + — > t + b + g is given by 

a s C F N c 



T(H + - tbg) 



4ir 2 mH+ 



(Vl + V2){ J 1 ~ ( m H+ ~ m t- m b) J 2 



(22) 



+ (m H+ - m 2 - m b ) I 12 } + km t m h yiy 2 {J2 - (m H + - m t - m b )I 12 }] , 



with the integrals 



'12 



Jl 



1 



4m 2 H+ L 



K 



■ 



+2Sp(l - 131) - Sp(l - ft) - Sp(l - Hi 

11 11 

2 1 2 2 2 1 2 2 



(ft + 6m t m b ) ln/3o — -K( /;.' 



/ 2 2 2\ 

(m^+ -m t -m b ) 



mpn + m\l 22 + Ji + h 

[2k In f Wrrnmt N + ^ + (<+ + ^ + ^ ^ 
/ tm H+ [ y k J 

+2m 2 t In ft + 2m 2 In /3 2 . 



(23) 



(24) 



(25) 



Here 



/3o 



771^+ — m 2 — m 2 + k 
2mtmb 
_ "T,ff+ + m 2 — m 2 — k 



Sp(a; 



— m 2 + m 2 — k 

2m H +m b 

x dt , ,„ 

-In , 
o i 



(26) 



and «; = K(m H+ ,m 2 ,ml). The definitions and the explicit forms of the J's are given 
in 
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The one-loop corrected decay width to 0(a s ) including the real gluon emission can 
be written as: 



r corr (# + -> tb + tbg) = T(H + -> *6) + r(# + -> %) 



r tree^+ _^ ^ + ^ r ( gluon ) + 5r(gluino) , (27) 
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with T^ ree given by eq. (Q), and 

5r(gluon) = ^3 [2{m 2 H+ - m 2 - m\) ( yi Re(6Y} 9) ) + y 2 Re(5Y 2 {9) )) (28) 

4m t m b (y 1 Re(5F 2 (9) ) + y 2 Re(5Y 1 {9) ))] + T(H + -> tig) , 



5r(gluino) = * c 1 \2{m 2 H+ - m\ - m 2 ) ( yi Re(5Yf 9) ) + y 2 Re(5Y 2 { ' 9) ) 

- Am t m b (yxRe(5Y 2 iS) ) + y 2 Re(6Y} S) ))] . (29) 

We have checked that the corrected width of eq. (|27|) is infrared finite. 

4 Numerical Results and Discussion 

We now turn to the numerical evaluation of the corrected width eq. (p7j) . As the 
standard QCD corrections have already been calculated J7|, it is interesting here 
to study the influence of the gluino (and ti,bj) exchange corrections <5r (gluino). 
The whole analysis depends on the following parameters defined at the weak scale: 
m H +, tan/5, /i, A t , A b , Mq, Mjy, M^, and m g . For simplicity we assume Mq = Mq- = 
Mp and A t = A b = A. We have found that our final results are rather insensitive 
to these assumptions. We take m t = 180 GeV, m b = 5 GeV, mw = 80 GeV, mz = 
91.2 GeV, sin 2 ^ = 0.23, #7(4vr) = «2 = a/sm 2 8 w = 0.0337 and a s = a s {m H+ ). 
We use a s (Q) = 12vr/{(33 - 2n f )\n(Q 2 /A 2 nf )} with a s (m z ) = 0.12 and the number 
of quark flavors rif = 5(6) for m b < Q < m t (for Q > m t ). 

In Fig. 2 we show the dependence of ST (gluino) as a function of A and Mq 
for tan/3 = 2 (a) and 12 (b), and (m H +,mg, //) = (400,550,300) (GeV). We see that 
the size of the SUSY-QCD correction <5r(gluino) can be large going up to ~ 50% 
and that it can be comparable to or even larger than the standard QCD correction 
<5r(gluon) in a large parameter region. For fixed tan/3, 5r(gluino) has a strong de- 
pendence on the parameters Mq and A which determine the masses and couplings 
of £i ; 2 and &i 2 - <5r(gluino) is smaller for larger masses of ti and bx- for tan/3 = 2, 
the correction due to 5r(gluino) is about -15% for (Mq,A) =(100 GeV, 300 GeV) 
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(where m~ tl ~ 119 GeV, and mg i ~ 98GeV), but it is still ~ —5% for larger squark 
masses (Mq,A) =(400 GeV, 300 GeV) (where m h ~ 405 GeV, and m~ hl ~ 399GeV). 
This tendency is consistent with the decoupling theorem for the MSSM. Notice also 
the different behaviour for tan/3 = 2 and tan/3 = 12. 

In Fig. 3 we show the m#+ dependence of r tree , r tree + 5r(gluon), and r corr = 
r tree + £r(gluon) + 5r(gluino) for tan/3 = 2 (a) and 12 (b), and (m s , //, Mq, A) = 
(400, —300, 200, 200) (GeV). The parameter values correspond to fixed stop and sbot- 
tom masses: mi x = 90 GeV, m t - 2 = 366 GeV, m bi = 193 GeV, and m~ b2 =213 GeV (for 
tan/3 = 2) and m t ~ = 173 GeV, m h = 333 GeV, m~ bl = 152 GeV, and mg 2 = 247 GeV 
(for tan/3 = 12). (Note that for m g = 400 GeV the DO mass limit of the mass- 
degenerate squarks of five flavors (excluding ti j2 ) is rrig > 140 GeV ]T4| .) We see 
again that the correction <5r (gluino) can be quite large and that it is comparable 
to or even larger than 5r(gluon) in a large region. Quite generally, the corrections 
<5r(gluon) and <5r(gluino) are bigger for larger tan/3, but it can happen that they 
partly cancel each other. The correction 5r(gluon) has already been calculated in 0. 
Our results on <5r(gluon) agree numerically with ref. 7 within 10%. 

In Fig. 4 we show a contour-plot for ^ r jjf c 1 < ^ 1 r no ' > in the tan f3 — m- g plane for 
(m#+, /i, Mq, A) = (400,-300,250,300) GeV. This correction rises with increasing 
tan/3, going up to 50%! Concerning the rrig dependence, ^'f 1 ™^ increases up to 
m g = 300 — 450 GeV and then decreases gradually as m g increases. It is striking that 
even for a large gluino mass (~ 1 TeV) ST( flv™^ is larger than 10% for tan/3 > 3. 
From Figs. 2 and 4 we see that the correction 5T (gluino) decreases much faster for 
increasing Mq than for increasing m g . 

In Fig. 5 we show contour lines of <5r(gluino) in the //— A plane for tan /3 = 2 (a) 
and 12 (b), and (m H +,m g , Mq) = (400,550,300) GeV. This correction has a strong 
dependence on /i and a significant dependence on A. We have found that the 
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A dependence for tan/3 = 1 is much stronger than that for tan/? = 2. 



The reason for the large contribution of <5r(gluino) as compared to 5r(gluon) 
is the following: The vertex-correction part of the gluino-exchange [gluon-exchange] 
corrections (see Fig. lb and eq. (5)) is proportional to the H + tb coupling [H + ib 
coupling] which is essentially ~ (A t + /4tan/3)/t t cos/3 + (At, + fi cot j3)hb sin/3 [~ 
h t cos (3 + hb sin /3] . Hence the vertex-correction part of the gluino-exchange cor- 
rections <5r(gluino) can be strongly enhanced relative to that of the gluon-exchange 
corrections <5r(gluon) in the case the g-mixing parameters A and fi are large. In this 
case ti and b\ tend to be light due to a large mass-splitting. Note that the 6-mixing 
effect plays a very important role for large tan f3. 

5 Conclusion 

Summarizing, we have performed a complete calculation of the 0(a s ) QCD correc- 
tions to the width of H + -> tb within the MSSM. We have found that the QCD 
corrections are quite important. A detailed numerical analysis has shown that the 
SUSY-QCD corrections (due to gluino, t and b exchanges) can be comparable to or 
even larger than the standard QCD corrections in a large region of the MSSM pa- 
rameter space; here the mixings oHi — tn and &l ~ bn play a crucial role. This could 
significantly affect the phenomenology of the H + search. 



After having finished this study, we have been informed on a recent paper []I5] dealing 
with the same subject. 
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Figure Captions 



Fig. 1 All diagrams relevant for the calculation of the 0(a s ) QCD corrections to the 
width of H + -> tb in the MSSM. 



Fig. 2 Contour lines of #r(gluino) (GeV) in the A — Mq plane for tan/5 = 2 (a) and 
12 (b), and (m H +,m g , /x) = (400,550,300) (GeV). For these parameter values 
one has (r tree (GeV), <5r(gluon) (GeV)) = (4.10,0.31) and (1.91,-0.66) for 
Figs. 2a and 2b, respectively. The shaded area is excluded by the LEP bounds 
m h 6 ~ 45 GeV. Note that for m g ~ 550 GeV one has no squark mass bound 
from DO experiment |T4 . 



Fig. 3 rriH+ dependence of r tree (dashed line), r tree + 5r(gluon) (dot-dashed line), and 
r corr = r tree + ^(gluon) + <5r(gluino) (solid line) for tan/? = 2 (a) and 12 (b), 
and (m~ g ,Li,MQ,A) = (400,-300,200,200) (GeV). 



Fig. 4 Contour lines of <5r(gluino)/T corr in the tan f3-m g plane for (m#+, //, Mq, A) = 
(400, -300, 250, 300) (GeV). The area below the dotted line is excluded by the 
LEP limit m-+ > 45 GeV (assuming rrig = (a s / 1 a<i)Mi ~ 3.56M 2 ), where 
a 2 = (; 2 / (47r), M 2 is the SU(2) gaugino mass, and m-+ is the lighter chargino 
mass. 
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5 Contour lines of 5r(gluino) (GeV) in the \x — A plane for tan/3 = 2 (a) and 
12 (b), and (m H +,rrig, Mq) = (400,550,300) GeV. For these parameter values 
one has (T tree (GeV), 5r(gluon) (GeV)) = (4.10,0.31) and (1.91,-0.66) for 
Figs. 5a and 5b, respectively. The shaded area is excluded by the LEP limits 
mr 7 -+ > 45 GeV. For m- a ~ 550 GeV one has no limit from the DO 



experiment [14 



12 




Fig. lc 




Fig. Id 
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